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PREFACE 


The  work  described  in  this  publication  was  sponsored 
by  the  United  States  Army  Intelligence  Center  and  School. 
The  writing  and  publication  of  this  paper  was  supported 
by  the  Jet  Propulsion  Laboratory,  California  Institute 
of  Technology,  under  a  contract  with  the  National 
Aeronautics  and  Space  Administration,  NAS  7-918,  RE  182, 


EXECUTIVE  SUMMARY 


This  Technical  Memorandum  was  prepared 
originally  as  part  of  the  Generic  ELINT/COMINT  Sensor 
Report  (FY-86)  which  was  elmininated  under  the  FY-87 
statement  of  work  (SOW  #2)  ,  undated  (delivered  to 
JPL  19  November  1986) . 

The  purpose  of  the  Generic  ELINT/COMINT  Sensor 
Report,  of  which  this  paper  was  intended  (in  its 
final  form)  to  become  part  of,  was  to  establish  a 
basic  superhetrodyne  receiver  based  sensor  model  and 
perform  simulations  with  it  to  determine  the  shaping 
or  coloring  of  the  statistical  distributions  of  the 
radar  free-space  signal  parametrics  by  a  typical 
sensor  prior  to  reaching  the  self-correlation 
processes.  It  was  also  intended  for  incorporation 
into  the  algorithm  test  bed  so  algorithms  could  be 
tested  with  realistic  distorted  data  rather  than 
unrealistic  stastically  pure  data. 

This  work  was  originated  in  support  of 
unanswered  questions  from  previous  self-correlation 
studies.  The  modeling  and  simulation  approach  was 
used  because  "live  date"  could  not  be  obtained. 

This  paper  is  being  published  because  it  was 
completed  in  FY-86  with  FY-86  funds  and  still  serves 
a  useful  function. 

Several  significant  results  were  noted  while 
performing  simulations  on  this  model.  The  first 
order  math  model  is  developed  for  the  Bragg-cell 
receiver.  Whether  the  input  distributions  were 
uniform  or  Gaussian,  the  output  distributions  were 
found  to  be  shaped  so  as  to  append  a  staircased 
triangular  distribution. 
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ABSTRACT 


Bragg-cell  receivers  are  employed  in  specialized  Electronic  Warfare 
(EW)  applications  for  the  measurement  of  frequency.  Bragg-cell  receiver 
characteristics  are  fuily  characterized  for  simple  RF  emitter  conditions, 
but  less  understood  for  complex  and  wideband  RF  emitter  signals.  This 
receiver  is  early  in  its  development  cycle  when  compared  to  the  IFM  receiver. 

Functional  mathematical  models  are  derived  and  presented  in  this 
report  for  the  Bragg-cell  receiver.  Theoretical  analysis  is  presented  and 
digital  computer  signal  processing  results  are  presented  for  the  Bragg-cell 
receiver.  Probability  density  function  analysis  are  performed  for  output 
frequency. 

Probability  density  function  distributions  are  observed  to  depart  from 
assumed  distributions  for  wideband  and  complex  RF  signals.  This  analysis 
is  significant  for  high  resolution  and  fine  grain  EW  Bragg-cell  receiver 
systems. 
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APPENDIX  A 


The  development  effort  is  for  theoretical  analysis,  software  algorithms 
implementations,  signal  processing  analysis,  and  test  analysis  for  probability 
density  function  characterization  of  a  Bragg-cell  receiver  system.  SigPro 
Systems  Inc.  is  performing  this  work  under  JPL  Contract  No.  957474,  Mod. 
No.  1. 

This  final  report  describes  the  completed  study  on  the  Bragg-cell 
receiver.  A  functional  mathematical  model  is  derived  for  the  receiver. 
This  math  model  is  implemented  into  digital  computer  software  and  signal 
processing  analyses  are  performed. 

Theoretical  analysis  is  developed  for  the  Bragg-cell  receiver.  Ideal 
and  practical  subsystem  implementations  are  presented  and  analyzed.  The 
key  ELINT  parameter  of  interest  is  frequency,  hence,  the  theoretical  analysis 
presents  the  basic  analysis  approach  for  frequency.  Probability  density 
function  analysis  is  presented  for  first-order  subsystem  model. 

The  theoretical  model  for  the  Bragg-cell  is  developed  into  a  digital 
computer  program.  All  computer  programs  are  listed  in  Appendix  A. 

Signal  processing  analysis  is  performed  using  the  derived  first-order 
model.  The  pdf  characteristics  of  the  Bragg-cell  receiver's  output  frequency 
are  determined  for  assumed  input  frequency  pdf  (constant  and  Gaussian) 


BRAGG  CELL  RECEIVER 


INTRODUCTION 

A  Bragg-cell  receiver  is  used  to  measure  frequencies  of  input  RF  radar 
pulses,  with  wide  frequency  coverage,  large  dynamic  ranges  and  simultaneous 
or  time  coincident  signals.  The  ideal  Bragg-cell  receiver  will  measure  the 
frequency  parameter  on  a  single  pulse  basis  with  no  distortion.  Practical 
Bragg-cell  receivers  will  distort  the  frequency  characteristic  of  a  radar  emitter 
because  the  ideal  Bragg-cell  receiver  characteristics  are  difficult  to  approximate 
in  real  life. 

The  radar's  frequency  parameter  can  be  characterized  using  a  random 
variable.  The  random  variable  may  be  quantified  by  the  probability  density 
function  (pdf),  mean  value,  variance  and  other  moments.  As  a  radar  signal 
passes  through  the  Bragg-cell  receiver  the  random  signal  description  is  modified 
or  distorted  by  the  receiver's  nonlinear  transfer  function  characteristic. 
Bragg-cell  receiver  characterization  of  RF  frequency  will  differ  from  actual 
radar  emitter  characterization  by  multifaceted  distortion  effects. 

This  section  presents  a  functional  description  of  the  Bragg-cell  receiver 
with  potential  distortion  sources  and  mechanisms  and  a  first-order  Bragg-cell 
receiver  model  is  developed.  This  Bragg-cell  receiver  model  is  used  to  analyze 
the  distortion  of  the  frequency's  pdf.  The  first-order  Bragg^-ceJSl  receiver  model 
is  limited  to  key  distortion  effects. 

IDEAL  BRAGG-CELL  RECEIVER  CHARACTERIZATION 

A  functional  block  diagram  of  a  Bragg-cell  receiver  Is  shown  in  Figure 
1.  The  basic  receiver  includes  an  RF-to-IF  downconverter,  Bragg-cell  with 
AO  transducer,  laser  source,  beam  expander  optics,  beam  focusing  optics, 
detector  array  and  processing  electronics. 


RP-TO-IP 


The  RF  input  signal  is  characterized  by 


=  X&)  cesf^-t  + 


Eq.  1 


where 


Pulsed  or  CW  signal, 
amplitude  modulation  function, 
fixed  angular  carrier  frequency, 
phase  deviation  function. 


The  RF-to-IF  downconverter  is  described  in  Reference  1.  The  ideal 


downconverter  is  composed  of  an  ideal  local  oscillator  which  is  a  perfect  CW 
signal,  a  perfect  product  mixer,  and  an  ideal  bandpass  filter.  The  ideal  IF 
signal  output  from  the  downconverter  is  described  quantitatively  by: 


UjCtl  52.  CO$  UJxpti  -v*  4*0] 


Eq.  2 


Equation  2  is  identical  to  Equation  1  with  the  fixed  angular  frequency  being 
changed  from  to  *  .  The  IF  center  frequency  is  chosen  to  match  the 
center  frequency  of  the  Bragg-cell  and  AO  transducer.  Ideally  the  EF  signal 
replicates  the  input  RF  signal.  The  RF-to-IF  downconverter  is  described  and 
characterized  in  Reference  1.  RF-to-EF  downconverter  analysis  will  not  be 
presented  in  this  report. 

The  IF  signal  drives  the  AO  transducer  which  launches  an  acoustic  signal 
through  the  Bragg-celL  Specifically,  the  IF  electronic  signal  is  converted 
into  an  acoustic  wave  (slow  wave  replication  of  the  IF  signal)  which  propagates 


through  the  optically  transparent  3ragy-celL  Througn  the  elasto-optic  effect, 
the  acoustic  wave  produces  a  spatial  modulation  of  the  refractive  index  in 
the  Bragg-cell  (Reference  2.) 

The  laser  source  is  expanded  by  the  beam  expander  optics  and  illuminates 
the  Bragg-cell.  As  the  laser  light  passes  through  the  Bragg-cell.  the  refractive 
index  variations  produced  by  the  acoustics  wave  (which  is  a  slow  form  replica 
of  the  IF  signal!  are  impressed  onto  the  optical  signal  as  a  spatial  phase 
modulation. 

Focusing  optics  are  used  to  equivalently  Fourier  transform  the  modulated 
optical  beam.  Resultant  optical  Fourier  transform  signal  is  focused  onto  the 
detector  array  or  frequency  focal  plane  (see  Figure  1.)  The  photodiode  detector 
array  is  used  to  capture  and  convert  the  optical  signal  to  an  electronic  signal 
representation.  The  detector  array  is  a  linear  photodiode  array  which  is 
functionally  used  to  detect  and  measure  light  intensity  versus  frequency  of 
the  transformed  optical  signal.  The  detector  array  is  electronically  read  by 
the  processing  electronics,  which  also  provides  follow-on  signal  processing 
functions. 

The  Bragg-cell  modulator  is  depicted  in  Figure  2.  The  incident  laser 
light  interacts  with  the  acoustic  signal  in  the  cell  to  produce  a  diffracted  light 
information  signal  along  with  undeflected  laser  light  signal.  The  diffracted 
laser  signal  is  also  frequency  shifted  by  the  acoustic  signal  frequency.  The 
incident  laser  light  hits  the  Bragg-cell  at  an  angie  9^  with  respect  to  the 
*  axis.  The  IF  electronic  signal  is  sent  to  the  AO  transducer  to  convert  the 
IF  signal  to  a  replica  acoustic  signal.  The  acoustic  wave  of  frequency  f , 
propagates  with  velocity  v#  along  the  x  -  axis  in  the  Bragg-celL 

The  diffracted  laser  light  signal  results  from  the  AO  interaction  effects 
between  the  incident  laser  beam  and  the  acoustic  information  bearing  signal 


in  the  Fragg-cell.  Maximum  diffracted  laser  light  energy  occurs  in  the  first 
order  beam  when  (Reference  2): 


and 


a  Vs 


Eq.  3. 


Eq.  4 


where 

=  Wavelength  of  laser  light 

4",  =  Frequency  of  acoustic  wave 

\/s  =  Velocity  of  acoustic  wave  in  Bragg-cell  along  thu 

x-axis. 

©j  =  Effective  angle  of  incident  laser  light. 

9<i  =  Effective  angle  of  first  order  diffracted  laser  light 

information  signal. 

Detailed  Bragg-cell  analysis  would  require  AO  cell  refractive  index 
variations  analysis  for  precise  quantification.  This  first  order  Bragg-cell  model 
does  not  require  this  detailed  analysis.  Functionally  the  Bragg-cell  model 
is  accurate  and  conceptually  clear  for  the  external  angles  8  j  and  ed- 

Bragg-cell  processors  are,  in  general,  configured  so  that  8j<<  0.1 
radian,  hence,  Equation  3  can  be  readily  and  accurately  simplified  to: 


of  the  first  order  diffracted  laser  ligtit  information  signal  is: 


2  Vs 


Eq.  6 


Equation  6  reveals  that  9d  varies  directly  with  the  acoustic  wave  frequency. 
As  stated  earlier,  the  EF  signal  frequency  is  related  to  the  acoustic  wave 


where  F  is  the  focal  length  of  the  focusing  optics. 


A  fixed  Bragg-cell  receiver  will  result  in  Equation  9  reducing  to: 


<±  =  K,  ^ 


Eq.  10 


The  actual  displacement,  d,  of  the  diffracted  laser  beam  relative  to 
the  undeflected  laser  light  passing  through  the  Bragg-cell  is  directly  proportional 
to  the  IF  frequency  of  the  input  IF  electronic  signal.  The  diffracted  laser 
beam  will  hit  the  linear  detector  array  at  a  specific  location  which  corresponds 
to  one  IF  frequency.  The  Bragg-cell  receiver  is  calibrated  by  varying  the  IF 
frequency  of  the  input  IF  signal  over  the  entire  operating  bandwidth  and 
capturing  and  analyzing  the  electrical  output  signal  from  the  linear  detector 
array.  Detector  array  element  signal  versus  IF  frequency  characteristic  is 
used  to  quantify  the  Bragg-cell  receiver's  frequency  measurement  performance. 

The  undeflected  laser  light  is  (the  zero-order  output  from  the  Bragg-cell) 
described  by: 


SLqL±)  =  utfc)  casjVv^  +•  orcoj 


Eq.  11. 


where 


undeflected  laser  light  signal. 

*  amplitude  or  intensity  of  undeflected  laser  light  signal, 
a  laser  light  frequency 
OrCO  *  phase  variations  on  laser  light  signal. 

The  diffracted  laser  information  signal  is  quantified  mathematically 


.*«  * «  -  »  *  .  • 
V  •  _•  ", 


E q.  12. 


L/t')  cos 


where 

I  (£)  =  Diffracted  laser  information  signal. 

=  Amplitude  or  intensity  of  diffracted  laser  information  signal. 

*  Frequency  of  information  signaL 
<fe(£)  *  Phase  variations  of  information  signaL 

Equation  12  reveals  that  the  diffracted  laser  output  signal  is  frequency 
shifted.  The  frequency  shift  is  produced  by  the  laser  light  signal  and  acoustic 
signal  interaction  in  the  Bragg-celL 

The  Bragg-cell  receiver  is  conceptually  a  parallel,  multi-channel  spectrum 
analyzer,  which  determines  the  Fourier  magnitude  spectrum  of  the  input  IF 
signaL  The  electronic  circuit  model  for  the  Bragg-cell  receiver  is  depicted 
in  Figure  3.  The  IF  signal  input  is  applied  to  a  parallel  narrowband  filter  bank 
with  an  associated  diode  detector  and  Alter  bank.  Each  output  signal  is  a 
power  or  energy  Indicator  of  the  IF  signal  spectral  content  in  each  selected 
narrowband  Alter. 


PRACTICAL  BRAGG-CELL  RECEIVERS 
The  ideal  Bragg-cell  receiver  is  characterized  in  the  previous  section, 
with  the  measured  frequency  being  derived  from  Equations  9,  10,  and  12. 
This  section  will  identify  and  briefly  present  some  key  differences  between 
practical  Bragg-cell  receivers  and  the  ideal  bragg-cell  receiver.  The  Bragg-cell 
receiver,  shown  in  Figure  1,  is  a  reference  for  this  discussion. 


Electronic  Electronic 


Figure  3.  Electronic  Circuit  Model  for  Bragg-cell  Receiver 


BRAGG  CELL 


The  ideal  Bragg-cell  is  assumed  to  have  an  undistorted  acoustic  signal 
replica  (at  acoustic  frequency)  of  the  IF  input  signal  propagating  through  the 
cell.  The  AO  transducer  provides  this  electronic  IF  signal  to  acoustic  signal 
conversion.  The  AO  transducer  has  transfer  function  nonlinearities,  which 
produces  unwanted  and  added  acoustic  signal  components.  Simple  sinusoidal 
electronic  signals  are  subject  to  minor  variations,  while  complex  and  wideband 
signals  may  be  significantly  distorted  by  the  nonlinear  transducer  action. 

The  ideal  Bragg-cell  is  assumed  to  provide  a  perfect  acoustic  signal 
termination,  at  the  far  end  of  the  cell.  In  general,  the  acoustic  signal 
termination  is  not  perfect,  which  cause  low  level  acoustic  signal  reflections 
and  distortion  generated  signal  components  to  propagate  back  through  the 
cell.  These  reflected  acoustic  signals  are  considered  as  additional  distortion 
sources,  hence,  variations  in  the  diffracted  laser  information  signal  can  be 
expected. 

Actual  distortion  effects  are  determined  by  AO  transducer  nonlinearities, 
reflected  acoustic  signal  levels,  total  acoustic  signal  parameter  descriptors 
and  Bragg-cell  properties.  The  complete  acoustic  signal's  instantaneous 
frequency  is  determined  by  amplitudes,  frequencies,  and  phases  of  desired 
acoustic  signal,  undesired  reflected  acoustic  signals  and  nonlinearly  generated 
and  unwanted  acoustic  signals.  Frequency  variations  produced  by  instantaneous 
frequency  changes  will  produce  angle  variations  in  the  diffracted  laser 
information  signal;  thus,  frequency  measurement  errors  in  the  Bragg-cell 
receiver. 


LASER  SOURCE 


The  laser  light  source  is  ideaily  assumed  to  be  constant  amplitude  and 
fixed  frequency  signal.  Actual  laser  sources  can  have  some  minor  frequency 
shifts.  Minor  frequency  shifts  will  results  in  small  angle  variations  in  the 
diffracted  laser  information  signaL  (See  Equation  3.) 

Frequency  variations  will  widen  the  effective  spotsize  and  change  the 
angle  of  the  diffracted  laser  beam,  which  means  the  light  energy  may  be  spread 
over  more  detector  elements  in  the  linear  detector  array. 

The  laser  light  signal  does  not  have  a  constant  intensity  or  amplitude. 
The  amplitude  is  quantified  by  deterministic  and  random  signal  components. 
These  amplitude  variations  will  resuit  in  diffracted  laser  information  signal 
intensity  level  variations  at  the  detector  plane. 

OPTICS 

Beam  expander  optics  and  diffraction  signal  focusing  optics  are 
fundamental  modules  in  the  Bragg-cell  receiver.  Beam  expander  optics  shape 
the  coherent  laser  beam  to  illuminate  one  entire  side  of  the  Bragg-celL  This 
beam  expansion  provides  a  long  path  interaction  between  the  acoustic  signal 
and  laser  light  within  the  Bragg-celL  Beam  expander  optics  and  laser  source 
are  positioned  so  that  the  laser  light  strikes  the  cell  at  the  Bragg  angle. 
Nonlinearities,  positional  variations,  and  rotational  variations  in  the  beam 
expander  optics  can  produce  amplitude  and  angle  variations  in  the  diffracted 
laser  information  signaL  These  errors  are  usually  very  small  for  well-designed 
optics. 

Focusing  optics  are  used  as  the  so-called  Fourier  transform  lens.  The 
Bragg-cell  Is  positioned  in  the  front  focal  plane  of  the  lens.  The  photodetector 


array  is  positioned  at  the  back  focal  plane  of  the  lense.  (Reference  3.)  Again, 
nonlinearities  and  location  variations  in  the  focusing  optics  variations  can 
produce  intensity  and  angle  variations  in  the  diffracted  laser  information  signal. 
These  errors  are  also  very  small  in  well-designed  optics. 

DETECTOR  ARRAY 

The  ideal  detector  array  is  assumed  to  be  a  long  linear  array  with  high 
frequency  resolution  and  sensitivity.  Practical  Bragg-ceil  receivers  use  long 
linear  photodiode  arrays,  where  typically  512  or  1024  photodiodes  are  closely 
spaced  along  the  array  length. 

Figure  4  presents  an  illustration  of  two  diffracted  laser  information 
signals  striking  the  photodetector  array  at  two  different  diodes  (two  different 
frequencies.)  In  actual  practice,  the  light  spotsize  for  each  beam  is  larger 
than  the  area  of  one  photodiode;  hence,  the  light  spills  onto  the  area  between 
detectors  and  also  onto  adjacent  detectors. 

Diffracted  laser  information  signal  angle  variations  (frequency  variations) 
produce  energy  spread  over  several  or  many  photodiode  cells.  Each  photodiode 
cell  is  used  to  determine  the  signal  energy  in  resolved  and  calibrated  frequency 
space.  The  photodetector  array  is  electronically  readout  to  determine  frequency 
and  amplitude  information  for  each  diffracted  laser  light  information  signal. 
The  photodetector  array  is  functionally  and  practically  described  as  a  frequency 
sampling  unit,  with  512  to  1024  frequency  bins.  Each  photodiode  determines 
the  elemental  frequency  resolution. 

Bragg-cell  receivers  can  accurately  measure  frequencies  of  time  coincident 
or  time  overlapped  signals,  which  is  a  big  advantage  over  other  receivers, 
such  as  the  IFM  receiver.  These  receivers  are  available  with  bandwidths  to 
approximately  l  GHz,  and  frequency  resolutions  of  100  KHz  to  10  MHz.  The 


3ragg-ceil  receiver's  dynamic  range  is  approximately  30  to  40  db.  .vith  the 
photodetector  array's  dynamic  range  being  the  critical  limiting  factor. 

PROBABILITY  DENSITY  FUNCTION  ANALYSIS 
BRAGG-CELL  RECEIVER 

The  prooaoility  density  function  (pdf)  of  the  Bragg-cell  receiver's  measured 
frequency  parameter  is  considered  in  this  section.  The  pdf  characteristic 
is  developed  for  a  first  order  Bragg-cell  receiver  model  under  high 
signal- to-noise  ratio  conditions.  The  IF  signal's  frequency  characteristic  is 
assumed  to  have  a  constant  pdf  or  a  Gaussian  pdf. 

Ideal  and  practical  Bragg-cell  receivers  are  discussed  earlier  in  this  report, 
and  characteristics  equations  are  provided  for  the  ideal  receiver.  The  Bragg-cell 
and  AO  transducer  are  key  producers  of  first-order  distortion  effects.  Distortion 
effects  include  AO  transducer  nonlinearities  and  reflected  acoustic  signals 
returning  from  the  far  end  of  the  Bragg-celL  This  analysis  assumed  that  acoustic 
signals  generated  by  AO  transducer  nonlinear: ties  are  significantly  larger  than 
reflected  acoustic  signals. 

The  ideal  Bragg-cell  receiver  contains  a  single  acoustic  signal 

'jet)  =■  -t- d.Ct-fj  Eq.  13. 

which  results  in  the  diffracted  laser  information  signal  given  in  Equation  12. 
fjis  the  frequency  of  the  acoustic  signal,  which  is  directly  related  to  the 
frequency  of  the  input  IF  signal.  (See  Equation  7.) 

AO  transducer  nonlinearities  are  producing  undesired  acoustic  signals 
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in  the  Bragg-cell.  Undesired  or  unwanted  acoustic  signals  are  described  as 
a  cosinusoidal  series  of  N  components  (harmonics  and  intermodulation 
components  at  a  significant  energy  level)  which  are  characterized  as: 


N  r  — 

u  tf)  =  2.  ^  Cos  zv  ^  +  '*'*  ^ 
m=.  2  L.  _ 


E  q.  14. 


where 

Uit)=  Total  unwanted  acoustic  signals  in  Bragg-celL.  The  AO 
transducer  will  generate  harmonics  of  and  other 

intermodulation  frequency  components. 

U,«  =  Amplitude  of  n  th.  unwanted  acoustic  signal  component. 

3  Frequency  of  n  th.  unwanted  acoustic  signal  component. 

Harmonics  of  f1  and  intermodulation  frequency  components 
are  present. 

C 3  Phase  variations  of  n  th.  unwanted  acoustic  signal  component. 
The  total  acoustic  signal  in  the  Bragg-cell  is  found  by  summing  Equations  13 


and  14: 


<juC-£)=  3  Ob')  + 


Eq.  15. 


Equation  15  can  be  rewritten,  with  a  lot  of  work,  as: 
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Eq.  16. 
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The  Instantaneous  frequency  variation,  resulting  from  3  ( t )  phase 
variations  In  Equation  16,  is 
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Eq.  17. 
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The  instantaneous  frequency  variation  is  a  very  complex  function  of 
amplitudes,  frequencies,  and  phases  of  all  unwanted  acoustic  signal  components 
in  the  Bragg-celL  Equation  17  can  be  rewritten  using  the  Fourier  series 
expansion  as  a  fixed  frequency  plus  a  series  of  harmonically  related  frequency 
components.  A  first-order  approximation  is  the  reduction  of  Equation  17  to 
a  fixed  frequency, 
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Eq.  18. 


f2  is  used  to  approximate  the  resultant  frequency  component,  which 
is  produced  by  unwanted  acoustic  signal  components  propagating  in  the 
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The  total  acoustic  signal  (Eq.  16)  can  be  approximated  as: 


Eq.  20. 


The  total  instantaneous  frequency  of  v(  t )  (Equation  20)  is: 


fy-t)  =  ?,  +  ^ 


Eq.  21. 


Equation  12  can  be  easily  changed  to  approximate  unwanted  acoustic 
signals.  The  modified  Equation  12  results  in  the  following  equation: 

-a,!*)  =  L>)  C0SJW&+4,  +$l')t  +  §'(t)] .  ^  22_ 

The  actual  displacement  of  the  diffracted  laser  information  signal  (Equation 
22)  relative  to  the  undeflected  laser  light  is 

,  F  ?-($,  +  *.') 

d  —  - — -  *  Eq.  23. 

V5 

Equation  23  reveals  that  d  is  a  random  variable  since  it  is  a  summation 
of  two  random  variables  f  1  and  f2.  Conceptually,  the  spot  size  of  the 
diffracted  laser  information  signal  is  increased  by  the  unwanted  acoustic  signal 
components.  The  actual  displacement  d  is  directly  calibrated  to  a  photodiode 


ceil  in  the  linear  detector  array,  which  corresponds  to  a  specified  frequency 


f  1  and  Z 2  are  assumed  to  be  inaependent  random  variables  for  this 
analysis.  Pdf  characteristics  of  Z  1  ard  z  2  are  assumed  to  be  either  constant 
or  Gaussian.  The  pdf  characteristic  of  d '  is  the  convolution  of  pdf  (f^) 
and  pdf  (f  ) ,  if  the  random  variables  are  independent  (Reference  4.) 

O®0 

p4?  (d‘)=  (dU.)  pif  (*,)  df,  Bq.  24. 


where 


J'a.  d  7 


Each  frequency  bin  in  the  linear  photodiode  array  exactly  corresponds 
to  ad  bin  coverage  in  displacement  space. 

The  output  frequency  of  the  detector  array  is 


?0=  K.  d 


Substituting  Equation  2S  into  Equation  23,  the  result  is: 


E q.  25. 
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The  Bragg-cell  receiver's  measured  output  frequency  can  be  stated  in 
terms  of  the  Input  IF  signal's  frequency  by  using  Equations  26  and  a  modified 
version  of  Equation  7.  The  output  frequency  is: 
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Eq.  27. 


=  The  equivalent  IF  frequency  of  unwanted  acoustic  signal 
components  approximated  as  acoustic  frequency  l  . 
tQ  is  a  linear  combination  of  two  random  variables  f  and  f^2 
,  which  are  assumed  to  be  independent.  The  pdf  (f_)  is  the  convolution 

w 

of  pdf  (f  )  and  pdf  (  2  ) ,  which  is  expressed  mathematically  as: 


(?o)  —  P^C^Zf)  ^Xf  • 


Eq.  28. 


The  pdf  ( f  Q )  characteristic  is  functionally  divided  into  sampled 
bins  by  the  photodiode  cells  in  the  linear  detector  array.  This  frequency  domain 
binning  process  is  equivalent  to  discrete  sampling  of  the  pdf  ( t  Q  ) 

characteristics  in  AfQ  bin  widths.  Af3  is  the  basic  frequency 
resolution  of  the  Bragg-cell  receiver  at  the  linear  detector  array.  The  actual 
pdf  characteristic  of  the  Bragg-cell  receiver's  frequency  parameter  is  a  discrete 
characteristic,  which  is  a  function  of  the  pdf  of  the  IF  signal  frequency  and 
pdf  characteristic  of  the  equivalent  unwanted  signal  components. 
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SIGNAL  PROCESSING  ANALYSIS 

INTRODUCTION 

Signal  processing  analysis  results  are  presented  in  this  section.  The 
Bragg-cell  receiver  model  is  implemented  on  a  Hewlett  Packard  Integral 
computer.  Computer  software  program  listings  are  contained  in  Appendix 
A  of  this  report. 

Pdf  analysis  of  the  Bragg-cell  receiver's  frequency  parameter  is 
performed,  assuming  constant  and  Gaussian  pdf  characteristics  for  the  input 
signal’s  IF  frequency.  Pdf  analysis  is  performed  for  selected  frequency 
resolutions  of  the  Bragg-cell  receiver.  Also,  statistical  analysis  are  performed 
on  output  pdf  signals  and  summary  results  are  presented. 

BRAGG-CELL  RECEIVER  MODEL  RESULTS 

The  Bragg-cell  receiver  is  math  modeled  in  the  previous  section. 
First-order  model  implementation  is  shown  to  be  a  convolution  of  pdf  (fjp) 
and  pdf  (f *)  for  the  Bragg-cell  receiver's  output  frequency  pdf  characteristic. 
Equation  28  is  the  final  descriptive  equation  for  the  output  frequency  pdf 
characteristic. 

Bragg-cell  receiver  signal  processing  results  are  summarized  in  Table 
1.  Signal  processing  analysis  programs  are  presented  in  Appendix  A.  Figures 
5  through  70  contain  detailed  pdf  plots  for  the  Bragg-cell  receiver's  output 
frequency  parameters.  These  plots  are  made  for  selected  pdf  characteristic 
of  f^  and  f  j  and  selected  Bragg-cell  receiver  frequency  resolutions. 

Table  1  defines  the  input  signal  or  fIF  in  terms  of  DF1  or  DF1  and 
The  pdf  (fjp)  is  assumed  to  be  a  constant  if  only  a  value  for  DF1  is  given. 
DF1  (MHz)  indicates  the  frequency  excursion  of  the  input  frequency  (f^) 

For  a  Gaussian  pdf  (fjp),  DFl(MHz)  and  o^(MHz)  is  selected  as 
1.166  MHz  or  4.166  MHz  for  Gaussian  pdfs. 
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The  distortion  signal  or  f  *  is  quantified  in  terms  of  DF2  or  DF2  and 
a,.  The  pdf(f^,)  is  assumed  to  be  a  constant  if  only  a  value  for  DF2  is  given. 
DF2(MHz)  indicates  the  frequency  excursion  of  the  distortion  signal  f^. 

For  a  Gaussian  pdf(f ^),  DF2(MHz)  and  c^MHz)  indicates  six  times 
sigma  and  standard  deviation,  respectively.  DF2  is  selected  as  2,  5,  9,  15, 
and  24  MHz.  c^MHz)  ^  selected  as  0.333,  0.833,  1.5,  2.5  and  4.0  MHz 
for  Gaussian  pdfs. 

Bragg-cell  receiver  frequency  resolution  (MHz)  is  chosen  to  be  0.5, 
1.0,  2.0  or  5.0  MHz.  These  frequency  resolution  selections  are  readily 
expected  with  current  Bragg-cell  receiver  technology. 

The  last  three  columns  of  Table  1  contain  statistical  results  of  the 
Bragg-cell  receiver's  output  frequency.  Sigma  (MHz),  skew  and  kurtosis 
are  computed  for  each  output  frequency  pdf  characteristic  shown  in  Figure 
5  through  70. 

The  pdfffj)  characteristic  is  a  staircased  trapezoidal  function  for  pdflf^p) 
=  constant  and  pdflf^)  =  constant.  Figures  5  through  22  show  numerous 
examples  of  pdf(fd)  characteristics  for  selected  input  signal  and  distortion 
signal  constant  pdf  characteristics.  Pdf(fd)  is  observed  to  significantly  depart 
from  a  constant  characteristic  and  approach  a  staircased  triangular 
characteristic  as  DF2  approaches  DPI. 

Figures  23  through  34  reveal  output  frequency  pdf(fd)  characteristics 
with  input  frequency  characteristics  assumed  to  be  Gaussian  pdf  and  distortion 
signal  frequency  characteristics  assumed  to  be  constant  pdf.  These  pdf 
characteristics  are  staircased  Gaussian  characteristics.  The  pdf  characteristic 
departs  from  the  Gaussian  characteristics  as  the  Bragg-cell  receiver's 
frequency  resolution  is  decreased. 


Figures  23  through  34  reveal  output  frequency  pdf(f^)  characteristics 
with  input  frequency  characteristics  assumed  to  be  Gaussian  pdf  and  distortion 
signal  frequency  characteristics  assumed  to  be  constant  pdf.  These  pdf 
characteristics  are  staircased  Gaussian  characteristics.  The  pdf  characteristic 
departs  from  the  Gaussian  characteristics  as  the  Bragg-cell  receiver's 
frequency  resolution  is  decreased. 

Figures  35  through  52  reveal  output  frequency  pdf(f^)  characteristics 
with  input  frequency  characteristics  assumed  to  be  constant  pdf  and  distortion 
signal  frequency  characteristics  assumed  to  be  Gaussian  pdf.  For  most 
example  plots  the  long  constant  pdf  characteristic  of  the  input  frequency 
tends  to  dominate  the  overall  pdf  characteristic.  The  pdf  characteristics 
are  approximately  a  staircased  trapezoidal  characteristic. 

The  pdfif^)  characteristic  is  a  staircased  Gaussian  function  for  pdf(f^p) 
=  Gaussian  and  pdf(f ~  Gaussian.  Figures  53  through  70  show  numerous 
examples  of  pdf(f^)  for  selected  input  signal  and  distortion  signal  frequency 
characteristics.  Pdf(fd)  is  observed  to  be  close  to  a  Gaussian  characteristic 
for  all  DF2  and  DF1  selections.  For  low  frequency  resolution  in  the 
Bragg-cell  receiver,  the  pdf  characteristic  departs  from  a  Gaussian 
characteristic. 
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CONCLUSIONS 


Bragg-cell  receivers  can  accurately  measure  frequencies  of  time 

coincident  and  time  overlapped  RF  signals,  which  is  a  big  advantage  over 

other  receiver  types.  Bragg-cell  receivers  are  available  with  bandwidths 
to  approximately  1  GHz,  and  frequency  resolutions  of  100  KHz  to  10  MHz. 

The  Bragg-cell  receiver's  frequency  parameter  is  analyzed  and  results 
are  presented  in  this  report.  Probability  density  function  analysis  and 

statistical  analysis  results  are  presented  for  this  receiver  for  selected 

frequency  resolution  capabilities. 

The  first  order  math  model  is  developed  for  the  Bragg-cell  receiver. 

The  input  signal's  frequency  is  assumed  to  be  a  constant  pdf  or  Gaussian 
pdf  characteristic.  The  primary  distortion  signal  is  generated  by  nonlinear 
acoustic  transducer  characteristics.  The  distortion  signal  is  also  assumed 
to  have  a  constant  pdf  or  Gaussian  pdf  characteristic. 

Staircased  pdf  characteristics  are  observed  for  all  pdf  frequency 
characteristics.  The  staircased  function  is  introduced  by  the  frequency 
sampling  of  the  photodiode  detector  array  used  in  the  Bragg-cell  receiver. 

The  output  frequency  pdf  characteristic  is  observed  to  depart  from  the  input 
signal's  frequency  pdf  characteristic  for  many  of  the  wideband  RF  signal 
cases  and  wideband  distortion  signals  produced  by  the  AO  transducer. 

Distortion  signals  with  large  DF2  values  can  significantly  change  the 

output  frequency's  pdf  characteristic.  Also,  the  frequency  sampling 

(or  detector  frequency  resolution)  at  the  detector  can  further  distort 

the  pdf  characteristic.  Low  frequency  resolution  at  the  detector  produces 

the  most  significant  distortions.  EW  emitter  classification  and  EW 

direction  finding  systems  can  be  affected  by  Bragg-cell  receiver  distortions. 

These  distortions  are  especially  important  for  wideband  radar  emitters. 
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APPENDIX  A 
BRAGG-CELL  RECEIVER 
CONVOLUTION  PROGRAM 


tO  !  ••BRAGS"  PROGRAM 

20  '  BRAGS- CELL  RECEIVER  ANALYSIS 

30  i  CONVOLUTION  OF  Ft  ARRAY  AND  FZ  ARRAY 

40  !  THIS  VERSION  PERMITS  STORAGE  OF  CONVOLVED  ARRAY  UITH  ITS  FREQUENCY  AXIS 

SO  !  STORAGE  TAXES  PLACE  PRIOR  TO  PLOTTING — A  STRING  NAMED  INFO*  IS  ALSO 

SO  !  STORED— SHOULD  CONTAIN  MINIMUM  INFORMATION  ON  STORED  DATA 

70  i  PLOT  30TH  INPUT  ARRAYS  ON  X  AXIS 

30  !  PLOT  CONVOLUTION  ON  FREQUENCY  AXIS 

30  !  LABEL  TOTAL  VALUES  FOUND  IN  CONVOLUTION  ARRAY 

100  !  ARRAYS  MAY  3E  ANY  COMBINATION  OF  FLAT/GAUSSIAN 

110  !  DECEMBER  IS8S 

120  i ISP  "DC  YOU  WANT  TO  SAVE  CONVOLVED  ARRAYS?"  3  LINPUT  "ENTER  Y/N“,XP* 

130  ! IF  KP*C1,1J*"Y“  THEN  30SU8  KEEP 

140  CIS?  "ENTER  INTEGER  LENGTH  OF  LONG  ARRAY"  3  INPUT  DF1 
ISO  !DI3P  "ENTER  LENGTH  OF  SHORT  ARRAY"  3  INPUT  DF2 
ISO  CIS?  "MAKE  THESE  ENTRIES  IN  MHZ" 

170  GI3P  "ENTER  DELTA  FREQUENCY  FOR  THE  LONG  ARRAY"  8  INPUT  OFF t 
130  013?  "ENTER  DELTA  FREQUENCY  FOR  THE  SHORT  ARRAY"  3  INPUT  DFF2 
130  GF2-GF1A' OFF l/CFFZ)  3  N-CFf* 2  3  3GF1 -GFF I  * 1000000/S  3  SGFL0-QFF2* 1000000/ S 
200  OISP  "OF  1  -" i GF 1 1  “  0F2-"  i  DF2  i "  N»"|N  3  NMR1  -1  /(  SGF1  *SQR<  2*PI ) ) 

210  OPTION  BASE  1  3  RAO  3  NMRLC«!/< SGFL0*SCR(2*PI> >  3  MU-500000000 

220  DIM  PXISOO) ,PH<  500) ,PY( 433 5 , TSt 323  ,T1*tS2J  ,FX(SOO) .XCSOO) ,Y(SOO> ,lNF0*t8O] 

230  LINPUT  "ENTER  «<80  CHARACTERS  DESCRIBING  DATA", INFO* 

240  REDIM  PYIN-1  >  ,FX<N-5  )  ,PX<N>  ,PH(N) 

2 SO  FOR  K»1  TO  N  !  FILL  INPUT  ARRAYS  WITH  ZEROES 
2 SO  PX<X>-0  3  PH<  X ) *0 
270  NEXT  K 

230  OISP  “  *****PICX  COMBINATION*-****" 

230  OISP  "  1 -BOTH  ARRAYS  FLAT" 

300  OISP  "  2 -SHORT  FLAT,  LONG  GAUSSIAN" 

310  OISP  "  3-LONG  FLAT,  SHORT  GAUSSIAN" 

320  OISP  "  4- BOTH  GAUSSIAN" 

330  OISP  “ENTER  YOUR  CHOICE"  8  INPUT  PIC 
340  IF  PIC-2  THEN  SHRT 
350  IF  PIC-4  THEN  GOSUB  AR2 
360  LSTEP-1 

370  FOR  X- 1  TO  N\Z  ! MAKES  FLAT  LONG  ARRAY 
380  PH(X)«l/(OFH) 

390  NEXT  K 

400  SHRT:  I  MAKES  FLAT  SHORT  ARRAY 

410  IF  PIC-3  THEN  GOSUB  AR3 
420  SSTEP-1 

430  FOR  X-N/4-CFZ/Z  TO  N/4+CFZ/2 
440  PX(K>»l/( DF2-1 ) 

4S0  NEXT  K 

460  IF  PIC-2  THEN  GOSUB  AR2 
470  CNVLV: 

480  YTOT-O 

490  FOR  K-1  TO  N- 1 


=00  YSUM-0 

SIO  STL *500- OFF  1  a  STPL*S30+0FFI  3  3T7L*<  STPL-STL>/<  N-1  )  a  IVAL-STL+ST2L/ 

523  FOR  J»1  TO  K 

530  YSUM«YSUM+PX<  J)»PH<X-J+1 ) 

540  PY< X >=YSUM/< 3TIL* 1000000) 

SSO  NEXT  J 

550  IF  X  MOO  30=0  THEN  OISP  •**<»•*  t K 
S70  NEXT  X 

530  FOR  X*  1  TO  N-1  !  MAKE  FREQ  ARRAY  FOR  XAXIS 
530  FX<  K  )*IVAL 
500  IVAL=IVALtSTIL. 

510  NEXT  X 
523  SUMY-0 

530  FOR  X= 1  TO  N-2  ! INTEGRATE 
540  XS*< FX(X  + I )-FX(X> >*1000000 
530  YS=<PY(X)+PYCX+1 > >/2 
530  SUMY-XS* YS+SUMY 
573  NEXT  X 

530  OISP  "SUM*" i SUMY 

530  I  IF  KP*tl.U*"Y"  THEN  30SUS  XP1 

730  CHOICE:  !  PICK  PLOT 

710  OISP  "  1-PLOT  BOTH  ORIGINAL  ARRAYS  SAME  AXIS’* 

723  OISP  "  2-PLOT  CONVOLUTION  ON  FREQUENCY  AXIS- 

733  OISP  "  3- EXIT  PROGRAM" 

740  OISP  -  4- STORE  OATA” 

7S0  OISP  “  ENTER  YOUR  CHOICE”  3  INPUT  CH 

750  IF  CH*4  THEN  G0SU8  XEEP 
770  IF  CH#1  THEN  730 

780  LINPUT  “EXTERNAL  PLOTTER?  Y/N“ .HP* 

730  IF  CH«3  THEN  BEEP  «  OISP  “DONE”  0  STOP 
300  IF  CH=1  THEN  G0SU8  PLT1 
813  IF  CH-2  THEN  GOSUB  PLT3 

823  PLOTTER  IS  1  Q  PEN  -1  8  GCLEAR  8  LOCATE  30.170.2S.80 
330  Y1=0  8  Y2*AMAXIY)*1 . 1  8  X1=STL  6  XZ-STPL 
340  START: 

350  SCALE  X1.X2.Y1.Y2 

850  FXO  2  8  LAXES  <  X2-X1  )/ 10.  <  Y2-YD/10.X1 .  Y1 
370  XAXIS  Y2.1X2-X1 )/10.X1 ,X2 
880  YAXIS  X2.CY2-Y1 >/10,Y1 ,Y2 

830  MORCAT:  1  ADO  MORE  OATA  SAME  SCALE 

300  IF  CH-1  THEN  350 
310  IF  CH»4  THEN  350 

920  DISP  "ENTER  BRASS- RECEIVER  FREQUENCY  RESOLUTION" 

930  INPUT  RES 

340  KRES-INT1 RES/I 2»0FF 1/N) ) 

950  FOR  K=1  TO  TOP 

950  PLOT  X( X ) ,Y( K )  8  IF  CH-1  THEN  1040 
970  IF  CH-4  THEN  1040 
980  FOR  L-1  TO  KRES 
930  IF  K+LWGP  THEN  1050 
1000  PLOT  XCX+D.YtK) 

1010  NEXT  L 
1020  X-K+L-2 
1030  L= 1 
1040  NEXT  K 

1050  IF  HP*-"Y"  THEN  EPLOT 
1060  IF  CH-1  THEN  GOSUB  PLT2 
1070  IF  CH-4  THEN  PL OUT 
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1383  IMAGE  “CONSTANT  POP  PGR  PI.  OF  1  *“  .00.3,“  CONSTANT  POF  FOR  FI  .  CFF2*“,CD.O 

1390  IMAGE  '‘CONSTANT  POF  FOR  FI.  CF1  .00.0.'*  GAUSSIAN  POF  FOR  FI  .  DFF2-“.D0.D 
•• 

113©  IMAGE  “GAUSSIAN  POF  FOR  FI.  OF  1 .00.0. “  CONSTANT  POF  FOR  FI  ,  CFFZ»“.OD. 
0.“ 

1113  IMAGE  “GAUSSIAN  POF  FOR  F 1 .  OF  1 »“ ,00.0, “  GAUSSIAN  POF  FOR  FI  .  0FF2-“.0D. 
0.'* 

1113  7S«“3RAGG-CSLL  RECEIVER  POF  ANALYSIS"  9  CSI2E  3.3,1 

1133  LORG  5  9  MOVE  XZ-<  X2-XU/Z ,  YZ+<  Y2-YD/9  9  LABEL  USING  “K"  i  T* 

1 14©  C3I2E  3 

11  S3  MOVE  X2-<  X2-X1 5/2 . Y2+<  Y2-Y1 )/ 1 3  9  IF  PIC-t  THEN  LABEL  USING  1383  5  OFF  1 .OFF 

1  9  C3I2E  3 

1 1  S3  IF  PIC-3  THEN  LABEL  USING  1393  i  OFF  1 .OFF!  9  "SIZE  4 

1173  IF  PIC-2  THEN  LABEL  USING  1133  \  0FF1 .OFFI  9  CSIZE  4 

1133  IF  PIC-4  THEN  LABEL  USING  1113  i  OFFI , OFF!  9  CSIZE  4 

1193  MOVE  X2-<  X2-X1 )/2 , Y1-<  Y2-Y1 )/5  S  LABEL  USING  “K“  »  “FO  IN  MHZ- 

1233  MOVE  Xt-<  X2-X1 )/7, Y2-<  Y2-Y1  )/Z  3  OES  9  LOIR  33  9  LABEL  USING  “K“  (  “Pdf<FD> 
X  1E3" 

1213  IMAGE  “TOTAL -“.00. 000 

1223  .MOVE  X*  ,  Y1  -<  Y2-Y1  )/S  9  LOIR  3  9  LORG  2  9  LABEL  USING  1213  j  SUMY 
1233  PLOUT; 

1243  LINPUT  “GUMP  GRAPHICS?  Y/N",OG» 

12=3  IF  CGS»“Y“  THEN  OUMP  GRAPHICS 
1233  EPLOT: 

1273  HP*»“N" 

1233  LINPUT  “PLOT  ON  EXTERNAL  DEVICE?  Y/N" ,PLE* 

1293  IF  PL£*»“Y"  THEN  PLOTTER  IS  70S  9  PEN  1  9  LOCATE  33.113,29.89  9  SOTO  START 
1333  IF  CH-4  OR  CH-2  THEN  CHOICE 
1313  BEE?  153.330  9  DISP  “OONE" 

1323  END 

1333  PLT1:  1  PLOT  SHORT  ARRAY 

1343  REDIM  XCN).Y(N) 

1353  IF  CH-1  THEN  MULT-1  ELSE  MULT-  ? 33333303 
1 3B3  FOR  K» 1  TO  N  I  FILL  X  ANO  Y  ARRAYS 
1373  X(K)«K  9  Y( X  )-PX( K I-MULT 
1383  NEXT  K 

1393  XI -1  9  XZ-N  9  Y1-3  9  YZ-AMAXI Y>*  1 . 1  9  TOP-N  9  LINE  TYPE  1 
1433  PLOTTER  IS  1  9  PEN  -1  9  SCLEAR  9  LOCATE  33,173,25,83 
1413  GOTO  START 
1423  RETURN 

1433  PLTZ:  1  PLOT  LONG  ARRAY  SAME  SCALE 

1443  T»-“ INPUT  ARRAYS" 

14S3  IMAGE  "CONSTANT  ARRAY  FOR  FI,  DF1-" ,00.0,"  CONSTANT  ARRAY  FOR  F2  ,  DFF2-" , 
00.00 

1483  IMAGE  "CONSTANT  ARRAY  FOR  FI,  OF1-“,DO.O,"  SAUSSIAN  ARRAY  FOR  F2  .  0FF2-" , 
00.00 

1473  IMAGE  "GAUSSIAN  ARRAY  FOR  FI,  OFT-" .00.0,"  CONSTANT  ARRAY  FOR  F2  ,  0FF2-" . 
00.00 

1483  IMAGE  "GAUSSIAN  ARRAY  FOR  FI.  OFT-", CO. 0."  SAUSSIAN  ARRAY  FOR  FZ  .  0FF2-" , 
00.00 

1493  LORG  5 

1S33  MOVE  XZ-<  XZ-XJ )/Z , Y2+<  YZ-Yt >/"3  9  IF  PIC-Z  THEN  LABEL  USIN6  1473  »  OFF  I. OFF 

2  9  CSIZE  4 

1513  MOVE  XZ-< X2-X1 )/2 ,Y2+< Y2-Y1 )/23  8  IF  PIC-I  THEN  LABEL  USING  1453  »  OFFI, OFF 
2  9  CSIZE  4 

1S23  MOVE  X2-<X2-XI )/2 , Y2+< Y2-Y1 )/"8  8  IF  PIC-3  THEN  LABEL  USING  t468  »  OFFI, OFF 
2  8  CSIZE  4 
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1  S3©  MOVE  X2-<  X2-X1  >/2 , Y2+<  Y2-Y1  )/Z0  3  IF  PIC-4  THEN  LABEL  USING  1430  ;  CFF1 ,QFF 

2  9  C3IZE  4 

1540  LORG  S  3  MOVE  XZ-<  XZ-X1  )/Z .  YZ+<  YZ-YD/3  9  LABEL  USING  “X"  i  T* 

1550  IF  CH-!  THEN  MULT*!  ELSE  MULT- i 000000000 

1SS0  MOVE  X1-1XZ-X1 )/7.Y2-<Y2-Y1 )/2  a  DES  9  LOIR  30  9  LABEL  USING  *X"  «  "VALUE  X 
" iMULT 

t  S70  FOR  X»t  TO  N 
1530  Y<  X  )-PH<  X  )-MULT 
1530  NEXT  X 

1S00  CH-4  a  LINE  TYPE  4  9  TOP-N 
1510  30T0  MOROAT 
1SZ0  RETURN 

1630  PL73:  !  FILL  CONVOLUTION  ARRAY 

1S40  REDIM  X(N-1) ,Y<N-1) 

1650  MULT* 1000000000 

1 650  FOR  K»t  TO  N-1 

1670  X( X  >«FX( X  5  9  Y< K )*PY<  X )*MULT 

1630  NEXT  X 

1 630  TOP*N- I 

1700  RETURN 

1710  ARZ:  ! SUB- GENERATE  GAUSSIAN  FOR  THE  LONG  ARRAY 

1720  Stt*flU- 3-SGF 1  9  3tp«rtU+3*SGF 1  9  STZ-< Stp-Stt )/(  N\2 )  0  IVL-Stt 

1730  FOR  X*l  TO  N\2 

1740  PH<  X 1-NMR1 *£XP( < -( IVL-MU ) *2 )/( 2-S6F1 *2  5)  9  PH< X )-PH<  X >-STZ 
1750  IVL-IVL+STZ 
1760  NEXT  X 

1770  IF  PIC-4  THEN  AR3 
1730  IF  PIC-2  THEN  CNVLV 
1730  RETURN 

1800  AR3:  1 SUB-SENERATE  GAUSSIAN  FOR  THE  SHORT  ARRAY 

1810  Stt-ftU- 3- SGFLO  9  Stp»f1U+3«SGFLQ  9  STZ-<  Stp-Stt )/CN/4+0FZ/Z-< N/4-0FZ/Z) )  9  I 
VL-Stt 

1820  FOR  X-N/4-QFZ/Z  TO  N/4+0FZ/Z 

1830  PX<  X )«NMRLQ*EXP< (t< IVL-MU)"Z )/( 2»SGFLQ“Z ) )  9  PX(  X J-PXC X )*STZ 
1840  IVL-IVL+STZ 
1850  NEXT  X 

1860  IF  PIC-3  OR  PIC-4  THEN  CNVLV 
1870  RETURN 

1880  XEEP:  I  FILENAME  ENTRY 

1830  01 SP  "FILENAME  MUST  BE  ENTERED  IN  UNIX  PATH  FORM  /OISCNAME/FILENAHE" 

1300  LINPUT  "ENTER  DESCRIPTIVE  FILENAME". F* 

1310  CREATE  FS. 1 .233*2+89 
1920  ASSIGNS  t  TO  F« 

1930  PRINT#  1  |  N.INFO«,FX( ) ,PY( ) 

1940  ASSIGN#  1  TO  • 

1950  SOTO  CHOICE 
1960  RETURN 


APPENDIX  A 
BRAGG-CELL  RECEIVER 

PDF  GENERATION  S  STATISTICAL  ANALYSIS  PROGRAM 


S  >8RAG3?0F  ’’HIS  PROGRAM  COMPUTES  3RAG3-CELL  RCVR.  PDF  ANO  STATISTICS 
; 3  :  -HIS  PRC3RAM  SETS  DATA  FROM  A  CONVOLUCLUTION  FILS-DOES  CEN  MOMS 

22  CATA  IS  "HEN  PLOTTED  UITH  LABELLED  CENTRAL  MOMENTS.  INT ,  ETC - 

33  OPTION  3ASE  : 

43  DIM  Fa<  S33 ; ,-H<  330  5 ,INFO*C  SOI .TSC7S3 

SC  CIS?  -PLEASE  USE  FILENAMES  UITH  0  OISITS  FOR  OF!" 

S3  LINPUT  -FILENAME?*  ,F* 

73  ASSIGNS  I  "0  ~S 
33  READS  t  s  N  3  N-N— ! 

30  REDIM  f\';n;,.=  H(N) 

•  03  REACS  :  i  INFO*  )  ,PH<  ' 


!  13 

-CR  X-» 

“A  »l 
«  W  N 

«  t* 

1  4 m*J 

FX :  :<  ;»fx 

C  X  '■*'  33 

i  33 

NEXT  K 

1  43 

AS3I3NS 

•  u  ♦ 

153 

CIS?  3  3IS?  INFO*  3  0I3P 

t  S3 

S-REU* 

(  r  *  /  9 

PS-siSC  1,13 

173 

IF  PS- 'A 

*  i-HEN 

PIO-I 

1  83 

IF  =S«-3 

*  "HEN 

PIC-0 

;ss 

IF  PS--C 

*  THEN 

PIC-3 

<■  UtJ 

IF  ?S-*j 

*  THEN 

PIC-4 

1 

L.  1  <J 

i  L  \ 

FSC  3 , 4 3 

;  9  dflg-ual<f*cs  .7  3 ; 

003  013?  "ENTER  RESOLUTION  IN  MHO*  3  INPUT  RES 
033  :<ST»INT(RE3/(0*QFI/N) ) 

043  0FM*AMAX<  FX !-AMlN<  Fa  5  9  FSTPM-QFM/< N-l  5*XST  9  CF-FSTPM* I  300030  9  SUMZ-0 
053  SUMY-3  3  SUMS-3  9  SUM3-0  9  SUM4-0 
033  FOR  K- !  TO  M-1 

273  XS-<FX?K+1  5-FX<K>>* 1000000  9  YS-< PH< K 5+PH( K+1 ) )/2 
233  3UM2»SUMZ+XS*YS 
233  NEXT  K 

333  FOR  K-l  TO  N  STEP  XST 

313  fx(k;»fx«k;* 1303300 

323  PH(X  )-PH<K 1/SUMZ 

330  SUMY-SUMY+PH<  K  )*FX.(  K  )*0F 

340  NEXT  K 

353  XMEAN-SUMY 

363  FOR  K»t  TO  N  STEP  XST 

373  3UMS-SUMS+<  FX 1 K 1-XMEAN  >  *2  »PH<  X  >  «QF 

380  NEXT  X 

330  UNC-SUMS  9  3I3MA-S0R< VNC > 

400  FOR  K-l  TO  N  STEP  XST 
410  SUM3-SUM3+<FX<K 1-XMEAN  )*3*PH< K  )*0F 
420  3UM4»SUM4+<  FX ( K  5-XMEAN  )*4*PH<  K  )  *0F 
433  NEXT  K 

440  CM3-SUM3  9  CM4-SUM4 

453  SKEU»CM3/VNC*1 .5  9  XURT-CM4,'UNC*Z 

483  PLOTTER  IS  1  9  PEN  -1  9  SCLEAR  9  LOCATE  30.170,25.30 

473  Y1 -3  9  Y2-AMAX(PH;*13Q00Q033*1  .1  9  X1-FXC 1 )/ 1000300  •  X2»XS«N/ 1 000000/2+530 
480  START: 


4 

*:u 

SCALE  ,:<2  ,yi  ,Y2 

533 

F'<0  3,:  3  LAXE3  1X2- 

<  1 

>  / 1 0 

Y2-YJ  :/i3 

,Y? 

■  ** 

2  1  *J 

X AX  IS  Y2 ,1  X2-<1  5/13, 

XI 

.XZ 

—  — 
Zm  4 

'•AXIS  <2  YZ-'I  5/13  , 

««  • 

:  i 

,Y2 

—  ^ 

MCROAT:  : 

AOO  MORE  OATA 

SAME  SCALE 

-  • 

2  **<; 

=CR  K— •  TO  N  STEP  X3T 

ZzZ 

IF  X+XSr-N  THEN  533 

zaZ 

PLjT  rXvK  /  /  i  4iu«/'u«0  ,r 

H( 

X  /*  t 

33300003  9  PLC 

1  P«lft  +  AJI  II  1  «w3luO0  , 

PHI  K 

;• 1 00003033 

Z  i  U 

=53 

NEX  T  K 

IMAGE  'CONSTANT  POF 

£ 

r  u 

R  FI 

,  OF!-*  ,300.0  , 

*  CONSTANT  P OF  FOR 

df=:-* ,000.0 

0 

5=3 

IMAGE  'CONSTANT  PDF 

fun  r  i 

r*C  •  -m  -  (A 

,  ur  1  *  ,  wUU  .  mi  , 

*  GAUSSIAN  POF  FOR 

r-  , 

0FF2-*  ,000.0 

mi 

3VU 

IMAGE  'GAUSSIAN  PDF 

^  0% 

r\j 

R  FI 

,  ur  I  *  ,m*uu  .  u  , 

*  CONSTANT  PDF  FOR 

"C  . 

off:-'  .odd. 

u 

313 

IMAGE  '3AU33IAN  POF 

FG 

R  FI 

,  0F1-*  .000.3, 

*  GAUSSIAN  POF  FOR 

“  i 

0FF2-'  ,000. 

529  T3-'3RAG6-CSLL  RECEIVES  PDF  ANALYSIS*  3  C3I2E  3.3,1 

533  LORG  5  3  MOVE  X2-<X2-:<1  5/2  ,Y2+<  Y2-Y1  )/ 9  3  LABEL  USING  *K*  t  T3 
**  **  —  ^  ^  • 

■a -t/  wiiiS  u 

5E3  MCVE  X2--X2-X1  5/2  ,Y2+<  Y2-Y 1  5/13  3  IF  PIC-1  THEN  LABEL  USING  =30  f  OF  I  ,DFL3  3 
*  C3IIE  a 

553  IF  PIC-3  THEN  LABEL  USING  530  «  DF1 ,3FL0  3  CSIZS  * 

5*3  IF  PIC-Z  THEN  LABEL  USING  BOO  I  OF!  ,DFL0  3  C3IZE  A 

330  IF  PIC-4  THEN  LABEL  USING  Slfl  t  0F1  ,OFLO  3  C3IZE  * 

530  MOVE  XZ-(XS-X1  !/C  , Y 1  - ( Y2-Y1  >/S  3  LORG  5  3  LABEL  USING  *K*  j  *Fd  IN  MHZ* 

TOO  LCRS  S  3  C3IZE  A  3  LORG  S 

713  MOVE  X1-<X2-X1 )/7,Y2-<Y2-Y1)/2  3  0E3  3  LOIR  30  3  LABEL  USIN6  *K*  «  *PdfCFd) 

X  1 E3* 

720  LOIR  3  3  RAC 
730  IMAGE  ’INTEGRAL-*  ,0.000 
7*3  IMAGE  '“MEAN-*  .OCDOO.OCO 
750  IMAGE  ‘SIGMA-* ,000.000 
750  IMAGE  * SKEW- *  ,300. 000 
773  IMAGE  ’KURTOSIS-* ,00.000 
780  CSIZE  2.S  3  LORG  2 

730  MOVE  X1+<X2-X1  5/2S ,Y2-< Y2-Y1  5/10  9  LABEL  USING  730  t  SUMZ 

300  LABEL  USING  743  I  X MEAN / 1 000000 

813  LABEL  USING  753  »  SIGMA/ 1 000000 

823  LABEL  USIN6  763  t  SKEW 

333  LABEL  USING  770  i  KURT 

840  LABEL  USING  *K*  i  ‘RESOLUTION-* |RE3 I*  MHZ’ 

853  LINPUT  ‘OUMP  GRAPHICS?  Y/N*  ,GOMS 
863  IF  GOM*C 1 ,1 >*Y*  THEN  OUMP  GRAPHICS 
873  IF  EPL*-‘Y*  THEN  900 

380  LINPUT  ‘PLOT  ON  EXTERNAL  OEVICE?  Y/N* ,EPLS 

890  IF  EPLS-*Y*  THEN  PLOTTER  IS  70S  9  PEN  I  8  LOCATE  30,110,29.39  9  GOTO  START 
303  SEEP  303,403  9  OISP  ‘DONE* 

9)3  ENQ 


